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1. INTRODUCTION 



Consider the nonlinear equation 

-u" + Q{x,u)u = k'^u, xeR, (1.1) 

where A; is a real parameter, the prime denotes the derivative with respect to the spatial 
variable x, and Q{x,u) has the form 

oo 

Q(a;,«)-^gn(a:)«^, (1.2) 

n=0 

with each Qnix) being real valued, bounded, measurable, supported in [0,6] for a fixed 
6 > 0, and the series 

( sup \qn{x)\]u^ (1.3) 

being entire in u. The assumption that the series in (1.3) is entire is equivalent to assuming 
that its radius of convergence is infinite, i.e. 

l/n 

Im I sup \qn{x)\ I = 0. (1.4) 



n— »-+oo 



a;€[0,6] 



Note that such an assumption is stronger than just assuming that Q{x,u) given in (1.2) 
is entire in u. 

In this paper we consider the solution to (1.1) satisfying 

u{0;k) = e, u'{0;k) = -ike. (1.5) 

As shown in Proposition 2.1, a unique solution to (1.1) exists for all a; e R and e R 

when l^l is sufficiently small. The results given in our paper hold both for real and complex 
values of e. We suppress the dependence of tt on £ for simplicity. Since for n > we assume 
Qn{x) = when x ^ [0, b], it follows that for each /c e R \ {0} the general solution to (1.1) 
for X ^ (0, b) is a linear combination of e*'^^ and e"*'^^. Thus, we have 

u{x;k) = se-''''=, x<0, (1.6) 
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which is equivalent to (1.5), and 

u{x;k) = A{k;e)e''''' + B{k;e)e-'^'', x>b. (1.7) 

Because of the resemblance with the (linear) Schrodinger equation, we will refer to 
(1.1) as a generalized nonlinear Schrodinger equation and to Q{x,u) as the (nonlinear) 
potential. The expressions for u{x;k) given in (1.6) and (1.7) indicate that we have a 
scattering problem in hand, where a plane wave is sent from x — +oo onto the nonhomo- 
geneity Q{x,u), and a part of the wave is transmitted to a; = — oo and a part is reflected 
back to a; = +oo. However, contrary to the linear case, there is in general no energy 
conservation, i.e. \B{k;s)\'^ — \A{k;e)\'^ in general depends on k. The linear part qo{x) of 
the potential Q{x,u) in (1.2) represents the restoring force density in wave propagation 
governed by (1.1) in the frequency domain. The higher order terms in the potential allow 
the description of nonlinear reaction of the medium to propagation of elastic waves. 

In analogy with the direct and inverse scattering problems for the linear Schrodinger 
equation, the corresponding problems for (1.1) can be formulated as follows. In the direct 
problem, given Q{x,u), our task is to determine the "scattering coefficients" A{k;e) and 
B{k; e) for sufficiently small |£|. On the other hand, in the inverse scattering problem, given 
some scattering data related to A{k; e) and B{k; e), the task is to recover Q{x,u). In our 
paper, we do not study the characterization of the scattering data so that the existence of 
a corresponding potential is assured; we only discuss the uniqueness and recovery aspects 
of our inverse problem by assuming that there exists at least one potential corresponding 
to our scattering coefficients. 

The inverse scattering problem analyzed in this paper is analogous to the recent study 
by Weder [1], where the time-dependent Schrodinger equation with a nonlinear term is 
investigated and both the linear and nonlinear parts of the potential are recovered from 
some appropriate scattering data by using a time-domain method. For other related studies 
of inverse problems on nonlinear equations using time-domain methods, we refer the reader 
to [2-5] and the references therein. 

Our paper is organized as follows. In Section 2 we show that the direct and inverse 

3 



scattering problems for the nonlinear equation (1.1) are equivalent to the corresponding 
problems for an infinite number of linear equations, and we analyze the basic properties 
of the scattering data for each linear equation and define the appropriate data set Vn 
given in (2.16) for each n > 1. We then solve the inverse scattering problem for each n > 1 
recursively. The solution of the inverse problem when n = 1 is well known, and in Section 3 
we list the basic facts from the case n = 1 that are needed later on to solve the inverse 
problems forn > 2. In Section 4 we show that the solution of the inverse scattering problem 
for each n > 2 can be obtained by inverting either of the two integral equations (4.4) and 
(4.5). In Section 5 we prove the unique invertibility of (4.4) and (4.5) to recover qn-i{x) 
for each n > 2, and we summarize the recovery of the nonlinear potential Q{x,u) in terms 
of the scattering data involving A{k; s) and B{k; e). Finally, in Section 6, we illustrate the 
direct and inverse problems for (1.1) with some concrete examples. 

2. PRELIMINARIES 

It is straightforward to verify that u{x; k) satisfies (1.1) and (1.5) if and only if it 
satisfies the integral equation 

1 

w(a;;/c) = £6"*'=^ + - / sin{k{x - t)) Q {t,uit;k)) u{t;k) dt. (2.1) 
k Jo 

Proposition 2.1. There exists a constant 6 > depending only on Q{x,u), but not on k 
(as long as k is real), such that if \e\ < S then a solution u{x;k) to (1.1) satisfying (1.5) 
exists for all a; e R, and it is unique. 

PROOF: Suppose that a solution u{x;k) to (2.1) ceases to exist, i.e. blows up, for some 
X e (0, b]. Then, for any r sufficiently large, there is an xq G (0, b] such that \u{xo; k)\ = r 
and \u{x; k)\ < r for x < xq. Fix one such r. Our assumption in (1.4) implies that there 
is a C > such that \Q{x, u) \ < C for x E [0, b] and \u\ < r. Let us take 

\s\<S:={r/2)e-^^\ (2.2) 
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Using I sin^l < |^| for real 9 and the realness of A;, for x < xq from (2.1) we get 

\u{x;k)\<\e\+ [ [x-t) \Q{t,u{t;k))\ \u{t;k)\dt 
Jo 

<\e\ + b / \Q{t,u{t;k))\ \u{t;k)\dt. 
Jo 

Applying Gronwall's inequality (see, e.g. Prob. 1 in Ch. 1 of [6]) on the last term above, 
we obtain 

\u{x;k)\ <\e\ + b\s\ J \Q{t,u{t;k))\ exp (^b Q{z,u{z;k)) dz^ dt. (2.3) 

Setting a; = a;o in (2.3), we have 

r<\e\ + b\e\ / C exp (Cb (xq - t)) dt < \e\ e^'' , 
Jo 

and a comparison with (2.2) indicates that 

< r < |£| e^^' < (r/2) e-^^'e^"' = r/2, 

which is impossible. Therefore, u{x; k) does not blow up in [0, b] and hence it exists for 
all a; e R. The uniqueness of u{x\ k) follows from the Lipschitz property of Q{x,u) with 
respect to u, which, in turn, follows from the analyticity of Q{x, ■). | 

From the conditions in (1.5), with the help of Theorem 8.4 in Sec. 1.8 of [6], when k 
is real and bounded, a; e R, and |£| is sufficiently small, we see that u is analytic in s and 
hence 

oo 

u{x; k) — ^ e^Un{x] k) — e wi(x; k) -\- U2{x-^k) + (2.4) 

n=l 

Note that the £°-term is absent in (2.4) because (1.5) implies that u{x; /c) = if £ = 0. 

We observe from (1.7) and (2.4) that the analyticity of m in £ for a; > 6 implies that 
A{k\ e) and B{k\ e) are analytic in £ at e = for real, nonzero, and bounded k. Thus, we 
have the expansions 

oo oo 

A{k- e)^Y. ^"^n(/c), B{k- e)^Y. ^^^^nik), (2.5) 

n=l n=l 



where we emphasize that the £°-terms are absent. With the help of (2.1) and the expan- 
sions (2.4) and (2.5), one can show that the generahzed nonhnear Schrodinger equation 
(1.1), the condition (1.6), and the expression (1.7) are equivalent to an infinite number of 
scattering problems for linear differential equations. The use of (2.1) allows us to avoid 
the interchange of the ^-differentiation and the infinite summation in (2.4). The resulting 
scattering problems for linear equations are, for n = 1, 

—u'l + qQ{x) ui = k^u\, a; e R, (2-6) 

wi(a;;/c) = e-^'^^, a; < 0, (2.7) 

ui{x;k) = Ai{k)e'^'' + Bi{k)e-~^^'', x>b, (2.8) 

and, for n > 2, 

-Un + qo{x)un ^ k'^Un - gn{x;k), xeR, (2.9) 

Unix;k) = 0, x<0, (2.10) 

«n(x;A;) = An(A;)e^'^" + S^(/c)e-^^^ x > b, (2.11) 
where we have defined 

gn{x;k) :^ qn-i{x)ui{x;k)'^ + hn{x;k), n>2, (2-12) 

n-l 

/i2(x; fc) := 0; hn{x; k) :^ '^Cjn{x; k) qj-i{x), n > 3, (2-13) 

with Cjn{x; k) for 2 < j < n — 1 being the coefficient of in the expansion of 

[eui{x; k) + e^U2{x; k) ^ h e'^~^Un-i{x; k)y . 

We list the first few hn{x; k) below: 

h2{x; k) = 0, h3{x; k) = 2uiU2qi, (2-14) 

h4^{x; k) — {2uiU3 + ul)qi + 'iu\u2q2, 
h^{x; k) = 2(ttitt4 + U2U3)qi + 2>{uiu\ + u\u3)q2 + Au\u2qz- 
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So far we have established the vahdity of Un{x;k) only for real and bounded k [cf. 
(2.4)] and those of A^ik) and Bn{k) [cf. (2.5)] only for real, nonzero and bounded k. In 
the next result, we analyze their analytic extensions in k to the entire complex plane C. 

Proposition 2.2. For n > 1, let Un{x;k), An{k), and Bn{k) he the quantities given in 
(2.4) and (2.5). Then, Un{x; k) for each x e R, kAn{k)., and kBn{k) are entire in k. 

PROOF: Since ui{x;k) satisfies the linear equation (2.6) with the condition in (2.7) and 
Un{x; k) for n > 2 satisfies the linear equation (2.9) with the condition in (2.10), it follows 
from Theorem 8.4 in Sec. 1.8 of [6] that Un{x; k) and tt^(a;; k) are entire in k. From (2.8) 
and (2.11), for a; > 6 and n > 1 we get 

Un{x] k) = An{k) e'^"" + Bn{k) e-^'^^, <(x; k) = ik An{k)e^^'' - ik Bn{k)e-''''' , 

and hence 

k Un {x; k) — iu'^{x; k) ^_ j^-p ^^^^ _ ^ '^n(x; k) -\- i u'^{x; k) ^^^^ (2 15) 

from which we see that kA^ik) and kB^ik) are entire in k. | 

Proposition 2.3. For each fixed k & C and all n > 1, the quantities Un{-;k) are bounded 
for X e [0, b]. For each fixed k & C and all n > 2, the quantities h„{- \ k) defined in (2.13) 
are bounded in x and vanish for x ^ [0, h]. 

PROOF: The boundedness of «„(•; k) for x e [0, b] can be established recursively for n > 1 
by using the fact that each tt„(-; A;) is a solution to a linear, ordinary differential equation 
with appropriate initial conditions at a; = [cf. (2.7) and (2.10)]. It is assumed that 
qn{x) for each n > is bounded in x and vanishes outside [0, h]. Since hn{x; k) is a linear 
combination of gi, . . . , qn-2 with coefficients that are polynomials inui, . . . , Un-i, it follows 
that hn{x; k) is bounded in x and vanishes when x ^ [0, 6]. | 

Recall that our aim in this paper is to solve the inverse scattering problem for (1.1), 
namely to recover Q{xj u) from some scattering data involving A{k; s) and B{k; e). In the 
light of (1.2) and (2.4)-(2.13), we see that our inverse problem is equivalent to the recovery 
of the qn-i{x) for each n>l from some data involving An{k) and Bn{k). We will establish 
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k A^i^k) 



the uniqueness of the recovery recursively; namely, first qo{x) will be shown to be uniquely 
recoverable from and then we will prove the unique recovery of Qn-i for 

each n > 2 from the data set V^, where we have defined 



As we will see, our uniqueness proof of recovery of Qn-i (x) will rely on the values of {k) 
and Bn{k) with large complex values of k. Recall that we have shown the validity of 
the solution u{x;k) to (1.1) with the condition in (1.6) only for real /c, and those of the 
scattering coefficients A{k;e) and B{k;e) only for real, nonzero, and bounded k. Hence, 
it is pleasantly surprising that we can prove the uniqueness of recovery of Q{x,u) from 
A{k;e) and B{k;e) without needing any extensions of the latter quantities to complex k 
values. 

3. INVERSE SCATTERING TO RECOVER qo{x) 

In order to solve the inverse scattering problem for (2.9), we need some basic facts 
related to (2.6). In this section we list those basic facts and refer the reader to [7-11] for 
details. Associated with (2.6), let L denote the unique selfadjoint realization of —cP/dx"^ + 
qQ{x) in L2(R). Even under weaker assumptions on qQ{x), namely, when qQ{x) is real 
valued, integrable, and vanishing outside [0,6], the following are known: 

(i) L has no positive or zero eigenvalues, it has no singular-continuous spectrum, and its 
absolutely continuous spectrum consists of [0,-|-oo). It has at most a finite number 
of (simple) negative eigenvalues, and we will denote the eigenvalues by — for j = 



(ii) The solution ui{x;k) to (2.6) satisfying (2.7) is usually known as the Jost solution 
from the right and sometimes denoted by /r(A;, x). As also indicated in Proposition 2.2, 
ui{x; k) is entire in k for each a; e R. From Lemma l(ii) on page 130 of [8], it follows 
that for /c G C+ and a; G [0, 6] we have 



X>„ := {A^(/c), B^{k), qo{x), 



qn-2{x)}, 



n>2. 



(2.16) 



1,...,7V. 




(3.1) 
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and for k e C+ \ {0} and x e [0, b] we have 



e^^"«i(a;;A;)-l| < ^ \qoit)\ dtj e^p (b \qoiz)\dzY (3.2) 



where we use C"*" for the upper half complex plane and put C+ := C"*" U R. Using 



(3.1) and (3.2), for k e C+ \ {0} and x e [0, b] we obtain 
|e^'=('^+i)^«i(x;/e)-l| < -^^^^ (^j^ \qo{t)\ d^j exp (^b{n + 1) \qo{z)\ dz^ . (3.3) 

(iii) Another solution to (2.6) which we denote by vi{x;k), usually known as the Jost 
solution from the left and sometimes denoted by f\{k,x), satisfies the asymptotic 
conditions 

wi(x;A;) = e^'=^[l + o(l)], v[{x; k) ^ ik e'''^[l + o{l)], x ^ +oo. 

Because of the compact support property of qo{x) we have 

r Bi{k) e'^'= - Ai{-k) e-^^'', a; < 0, 
v^{x;k)=l (3.4) 
[ e*'^^, x > b. 

For each a; e R, ui{x;k) is entire in k. The inequalities given in (3.1)-(3.3) also hold 
if we replace e^^^ui{x; k) in them by e~'^^^vi{x; k), which is a consequence of Lemma 
l(i) on page 130 of [8]. 

(iv) The scattering coefficients Ai{k) and Bi{k) given in (2.8) are related to the transmis- 
sion coefficient T{k) and the right reflection coefficient R{k) as Ai{k) = R{k)/T{k) 
and Bi{k) = 1/T{k). It is known that kAi{k) and kBi{k) are entire (of. Proposi- 



tion 2.2). The quantity Bi{k) is nonzero in C+ \ {mjjjLi and has simple zeros at 
k = iKj ior j = 1, . . . , N. In the so-called generic case Bi{k) has a simple pole at 
k = 0, and in the so-called exceptional case Bi{k) is continuous at A; = 0. 

(v) The potential qo^x) can be recovered from Ai[k) / B\[k) known for /c G R by any of 
the several methods [7-11]. No bound state data needs to be supplied [12-17] due to 
the compact support property of qQ{x), and in fact the bound state data (i.e. the 
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N constants Kj and the related norming constants) can be recovered via the unique 
meromorphic continuation of Ai{k) / Bi{k) from e R to e 0+ and using its poles 
and residues. Actually, knowledge of Ai{k)/ Bi{k) in some interval on the real /c-axis is 
sufficient to recover qo{x) because of the uniqueness of the meromorphic continuation. 

(vi) It is known that Bi{k) alone cannot uniquely determine qo{x), and unless B\{k) = 1 
there are infinitely many potentials corresponding to it. In the absence of bound 
states, i.e. if = 0, the coefficient Ai{k) uniquely determines qo{x) provided 1/Bi{k) 
vanishes at A; = 0. Otherwise, there are exactly two distinct potentials corresponding 
to Ai{k). When N > 1, there is no uniqueness and there are a discrete number of 
potentials corresponding to Ai{k). For details we refer the reader to [18]. 

(vii) As long as k'^ is not in the L2(R)-spectrum of L, i.e. when k e C+ \ {iK,j}jLi, the 
operator (L — A;^)"-^ is a bounded one from L2(R) to L2(R). The Green's function 
associated with (2.6), denoted by G{x,t;k), is defined as the integral kernel of the 
resolvent [L — k"^)^^ in the sense that for any g G iv2(R) we have 

/oo 
G{x,t;k)g{t)dt. (3.5) 
-oo 

It follows that, for any k e C"*" \ {iKj}jLi [cf. (i) and (iv) above], we get 

-Gxa:{x, t; k) + qo{x) G{x, t; k) = k^G{x, t; k) + 5{x - t), (3.6) 
where S{-) is the Dirac delta distribution. With the help of the Wronskian identity 

2ikBi{k) — v'i{x; k) ui{x\ k) — vi{x\ k) u\{x\ k), 
it can be verified that 



G{x,t; k) = 



2ik Bx{k) 
1 



vi{t\k)ui{x\k), x<t, 



(3.7) 



vi{x; k) ui{t; k), x > t. 



2ikBi{k) 

It is seen from (2.7), the second line in (3.4), and (3.7) that, when k e \ {iKj}jLi 
the Green's function G{x,t;k) decays exponentially as a; — > ±oo, and it is the only 
solution to (3.6) having this property. 
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4. INVERSE SCATTERING FOR qn-i{x) WITH n>2 



In this section we show that the recovery of qn-i{x) for each n > 2 is equivalent 
to inverting the integral equation given in (4.4) below or the one given in (4.5). Having 
seen in the previous section that qo{x) can be recovered uniquely from {Ai{k), Bi{k)}, we 
proceed recursively and prove the unique recovery of qn-i{x) from the data set Vn defined 
in (2.16). 



With the help of (2.6)-(2.11), we define 



yn{x; k) := Un{x; k) - ui{x; k), 

±>l{K) 



n>2. 



(4.1) 



The relevant properties of Vnix; k) are analyzed next. 

Proposition 4.1. Let yn{x; k) be the quantity defined in (4--1)- Then: 

(i) For each k e \ {iKj}^^]^, where {m-,}^^ is the set of zeros of Bi{k) in , the 
quantity yn{x; k) is a solution to (2.9), and it is the only solution belonging to L2(R) 
in X. In fact, yn{-', k) decays exponentially as x ^ ±oo. 



(ii) yn{x; k) satisfies 



yn{x; k) = < 



B^{k) 
' B^{k) 



-ika 



x < 0, 



(4.2) 



X > b. 



(iii) In terms of the nonhomogeneous term gn{x;k) defined in (2.12) and the Green's func- 
tion G{x,t;k) given in (3.7), we have 



/oo 
G{x, t; k) gn{t; k) dt, k e C+ \ {iKj}f^^. 
-oo 



(4.3) 



PROOF: First, from Propositions 2.2 and 2.3 and the fact that the only zeros of Bi{k) 
in C"*" occur at A; = inj for j — 1, . . . ,N, it follows that yn{x', k) is well defined for each 
k e 0+ \ {iK,j}jL-i^. Note that yn{x;k) solves (2.9) because Un{x;k) is a solution to the 
same equation and ui{x; k) is a solution to the corresponding homogeneous equation. We 
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obtain (4.2) directly from (2.7), (2.8), (2.10), and (2.11). For each fixed k e C+\{mj}^i, 
since e"*'^^ and e*'^^ decay exponentially as a; ^ — oo and as a; ^ +oo, respectively, we 
get the L2(R)-property of yn{-;k) stated in (i). Note that yn{x;k) is the only L2(R)- 
solution to (2.9) because the difference with any other L2(R) -solution must satisfy (2.6); 
however, as indicated in (i) of Section 3, the only L2 (R)-solutions to (2.6) occur when 
k = iKj for J = 1, . . .,N. Thus, we have proved (i) and (ii). From the properties listed 
in Proposition 2.3, it follows that gn{-;k) is bounded and supported in [0,6] and hence 
belongs to L2(R). In terms of the operator L defined in Section 3, yn{x;k) satisfies 
(L — k^)yn = —Qni and hence with the help of (3.5) we obtain (4.3). | 

In the next theorem we present the main integral equations from which qn-i{x) with 
n > 2 will be recovered. 

Theorem 4.2. For each n>2, the potential qn-i{x) satisfies 

f vi{t;k)ui{t;kyqn-i{t)dt = En{k), k e C, (4.4) 

/ m{t;kY+^qn-i{t)dt = Fn{k), keC, (4.5) 
Jo 

where En{k) and Fn{k) are completely determined by the data T>n defined in (2.16), and 
they are given as 

En{k) := -2ik Bn{k) - vi{t;k) hn{t]k) dt, (4.6) 

Jo 

F^{k):=2ik[B^{k)An{k)-A^{k)Br,{k)]- I ui{t;k) hn{t;k) dt. (4.7) 

Jo 



PROOF: Using (2.7) and (3.4) in (3.7), we get 
G{x,t;k) = < 



^ e-'''''vi{t;k), a;<min{0,O, 



2ikBi{k) 



{ 2ikBi{k) 

Using (4.8) on the right hand side of (4.3), we obtain 



e^'"'ui(t;k), x>max{b,t}. 



(4.8) 



yn{x;k) = 



rb 

„—ikx 



2ikBi{k) Jo 



vi(t;k) gnit;k) dt, x<0, 



^ "^'^^ / ui{t;k)gn{t;k)dt, x > b. 
Jo 



(4.9) 



2ik Bi{k) 
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Comparing (4.9) with (4.2), we see that 

2ikBn{k) = - f vi {t; k) gn{t; k) dt, n > 2, (4.10) 
Jo 

2ik[Bi{k) An{k) - Ai{k) Bn{k)] = f ui{t;k) gn{t;k) dt, n>2. (4.11) 

Jo 

With the help of (2.12), we rewrite (4.10) and (4.11) as (4.4) and (4.5), respectively. Even 
though we have derived (4.4) and (4.5) only for k G C"*" \ {iKj}jLi, with the help analytic 
extensions indicated in Proposition 2.2, we see that the former is actually valid for /c G C 
and the latter for A; G C \ {0}. Furthermore, the analytic extension to A; = for (4.5) is 
proved by showing that the quantity k[Bi{k) An{k) — Ai{k) Bn{k)] appearing in (4.7) is 
continuous at = 0. In order to see this, with the help of (2.15), as — > in C we obtain 

2ik[B,{k) Ar,{k) - A^{k) B^{k)] = <(x; k)[A,{k) + B,{k)] + 0(1). 

Although each of Ai (k) and Bi{k) may have a simple pole at A; = 0, their sum Ai (k) +Bi (k) 
is continuous aX k — because (2.8) implies 

Ai (k) + Bi (k) = e-^^'^m (x; k) + k B^ (k) ^ ~ ^, ^'^"^ , x>b, 

k 

and the right hand side isO(l)asA;^0inC. | 

With the help of Theorem 4.2, we see that our inverse scattering problem of recovery 
of qn-i{x) for each n > 2 can be stated as follows: Given E^ik) or Fn{k), determine 
qn-i{x) for X G [0, b] by inverting either (4.4) or (4.5). 

From (2.13) and (4.6) it is seen that the data set {i?„(/c), qo{x), . . . , qn-i{x)} uniquely 
determines all the quantities given in (4.4) except for q^-iix) there, and hence one can 
use (4.4) to determine qn-i{x) without needing An{k) in the data. 



5. UNIQUE RECOVERY OF qn-i{x) WITH n>2 



In this section we show that qn-i{x) with n>2 can be recovered by inverting either 
(4.4) or (4.5). We will present the inversion only for (4.5) because the inversion of (4.4) is 
similar. From (3.3) we see that we can write 



ui{t; k) 



n+l 



-ifc(n+l)t 



1 + 



U{t;k) 
k 



xG[0,6], A;gC+\{0}, 



(5.1) 
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where U{t; k) is uniformly bounded for x e [0, h] and k e C+; i.e. 



\U{x;k)\<M, a;e[0,6], /c e C+, (5.2) 

where M is independent of x and k. Let us separate the real and imaginary parts of the 
complex k variable in C"*", discretize the former, and keep the latter as a fixed positive 
parameter by writing 

^rn = I ^^"^s, + Z. (5.3) 

[n + 1)0 

Let us evaluate (4.5) at k — km and write the resulting equation in the operator form 
as K(f) = p, or equivalently as 

f K{m,t)(j){t)dt = p{m), m e Z, (5.4) 
Jo 

where we have defined [cf. (5.1) and (5.3)] 



K{m,t) := e-^^'^+i)*^! {t; kf^^ = e'^'^^"^*/'' 



\ ^ U{t; km) 



(5.5) 



27rm(n+ l)-i6-i + 

0(t):=e«(-+i)*g„_i(t), (5.6) 

p{m) -.^ Fn{km). (5.7) 

Note that we suppress the dependence on n and ^ because the inversion of (4.5) will be 
done at fixed values of n and ^. 

When qo{x) = 0, we have tti(t; k) = e"*'^* and U{t; k) = 0; in that case the operator 
K reduces to Kq, which is given by [cf. (5.5)] 

/■6 nb 

{Kog){m)^ / Ko{m,t) (l>{t) dt = / e'^'"''^*/^ (f>{t) dt. (5.8) 
Jo Jo 

The properties of K and Kq are analyzed next. 

Theorem 5.1. The operators K and Kq defined in (5.4) and (5.8), respectively, satisfy 

the following: 

(i) Kq maps 1/2(0, b) into ^2(2), and its operator norm is given by \\Ko\ \ = Vb. 
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(ii) Kq ^ exists as a map from ^2(2) into L2(0,6), and its operator norm is given by 
\\K^^\\ = 1/Vb. 

(iii) K — Kq maps L2{0,b) into ^2(2). The operator norm of K — Kq satisfies ||i^ — i^'o|| < 
•\/s(^), where s(^) is a monotone decreasing function on ^ E (0, +00) vanishing at 
infinity. 

(iv) There exists ^q > 0, determined by n, b, and \qQ{t) \ dt alone, such that the inverse 
operator K~^ exists for any C > Co- 

PROOF: For the proof of (i), note that Kq is the standard Fourier series operator and 
hence it maps L2(0, b) into ^2(2). With the help of the completeness relation 

^2M/h g-27rimx/6 = l).S{t-x), t, a; G R, 

mGZ 

using the definition of the operator norm, from (5.8) we obtain \\Kq\\ = \/b. As for (ii), 
the inverse Fourier series operator Kq^ is given by 

{K^'h){t) = \ J2h{m)e'-'-^'/\ (5.9) 

meZ 

and its operator norm is computed in a straightforward manner with the help of 

^2M/b ^-2^ijt/b dt = b- Smj, m, j e Z, 

'0 

where Smj denotes the Kronecker delta, yielding = 1/Vb. As for (iii), note that the 

square of the operator norm of K — Kq is defined as 



^0 



\K -K 



o| 



sup \\{K-KQ)gf= sup V|[(K-Ko)^]( 

lfll=i M-^mez 



m) 



where \ \g\ \ denotes the norm of g in L2(0, b). With the help of (5.3)- (5.5) and (5.8) we get 



\K — kq\\^ — sup 



mez 



/' 

^0 



U{t- km) 



2Trm{n + l)-^b-^ + i$, 



- (j){t) dt 



(5.10) 



Using the Cauchy-Schwarz inequality on (5.10), \\g\\ = 1, and (5.2), we obtain 



Uit; K 



27rm(n + l)-i6-i + 
1 



\K-KQf<Y, 



dt 



mGZ 
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and since the summation in the last term can be evaluated in a closed form, we obtain 

||;,_;,„||^<^l!q£±i)eoth((!i±l«), ,>0, n>2. (5.11) 

Note that coth(a;) is monotone decreasing on a; e (0,+oo). Using (3.3), (5.2), and (5.11), 
we see that ||K — KqH < •\/s(C)) where 

^(0 (^J^\qo{x)\dx^ coth(^^^^±^^ e^p (^2b{n + 1) Mx)] dx^ 

Next, we need to prove (iv). Writing 

K = Ko + {K-Ko)^ Ko[I + K-\K - Kq)], 
we see that can be evaluated as a Neumann series 

oo 

whose convergence is assured by choosing ||(i^ — Ko)|| < 1- Prom (ii) and (iii), we 

see that this is achieved by choosing as the unique solution to s(Co) = b. Then, for any 
^ > the existence of is assured. | 

Having shown that K is invertible for suflBciently large ^ values, we can recover qn-i{x) 
with the help of (5.4), (5.6), and (5.7) as 

qn-i{x) = e-^(^+^)^ {K-^Fn){x), x e [0,6], n > 2. (5.12) 

Let us now summarize the reconstruction in the solution of the inverse scattering 
problem for (1.1). Assume that we are given the scattering data consisting of A{k;s) and 
B{k;e) for all k in some subinterval of the positive /c-axis and for all |£| < 5, where 5 is 
some positive number (no matter how small). Our aim is to recover the potential Q{x,u) 
corresponding to the given scattering data. 

(i) Using (2.5), obtain An{k) and Bn{k) for n > 1. Note that kAn{k) and kBn{k) have 
entire extensions to C. 
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(ii) Recover qo{x) and ui{x;k) from {Ai{k), Bi{k)} by using any one of the available 
inversion methods described in (v) of Section 3. 

(iii) Recover qn-i{x) and Un{x;k) for n > 2 by using the data Vn given in (2.16) in 
a recursive way. The recovery of qn-i{x) is achieved by inverting either (4.4) or 
(4.5). The recovery of Un{x; k) amounts to solving the linear equation (2.9) with the 
condition in (2.10). 

(iv) Having recovered all the qn-i{x) for n > 1, we obtain Q{x,u) via (1.2). 

The uniqueness of the recovery of Q{x, u) from the data {A{k; e), B{k\ e)} is summa- 
rized next. 

Theorem 5.2. Consider the data consisting of A{k; e) and B{k; e) given for all k in some 
subinterval of the positive k-axis and for all \e\ < 5, where 6 is some positive number (no 
matter how small). Further assume that there exists a potential Q{x,u) corresponding to 
this data. Then Q{x,u) is the only potential corresponding to that data. 

6. EXAMPLES 

In this section we illustrate the direct and inverse problems for (1.1) in a special case 
with some explicit examples. Assume that Q{x,u) given in (1.2) has the form Q{x,u) = 
q2{x)u'^, i.e. qj{x) = for j > except when j = 2. From (2.6)-(2.8) we get Ui{x;k) = 
^-ikx xeK, Ai{k) = 0, Bi{k) = 1. From (2.9)-(2.13) with n = 2, we obtain U2{x; k) = 
for a; e R, A2{k) = B2{k) = 0. Using (2.9)-(2.11) with n = 3, we see that U3{x;k) 



with the initial conditions tt3(0; k) = u'^{0; k) = 0. Using variation of parameters, we get 



satisfies 



u'^(x; k) + k^us{x; k) = q2{x) e 



—Zikx 



x e R 





X 



U3{x;k) 



Q2{t) e 



-Aikt 



dt 



Q2{t)e 



-2ikt 



dt, 



xeK. (6.1) 



2ik 



Comparing (6.1) with (2.11) we see that 




(6.2) 
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Thus, we see that A^{k) and B^{k) are related to each other as A^{k) = —2B^{2k). With 
the help of e^^'^dk = 2T:5{a), the potential q2{x) is uniquely recovered from A^^^k) or 
Bz{k) given in (6.2) as 

g2(x) = -/ kA^{k)e''^''dk^ -- kB^{k)e^'^''dk. (6.3) 

Alternatively, ^'2(0;) can be recovered by using the method of Section 5 as follows. Let 
us suppose that we are given ^3 (A;). Using (2.14) and (4.7) we see that Fz{k) = 2ikA^{k). 
Thus, with the help of (5.3) and (5.7) we get 

p{m) = {irmv/b — 2^) A2{i^ + m7r/26), m e Z. 

Since qo{x) = 0, we have K = Kq, where K and Kq are the operators appearing in (5.4) 
and (5.8), respectively. Thus, for any ^ > 0, with the help of (5.9) and (5.12), we construct 
q2{x) explicitly on [0,6] via 

q2{x) = ^ e-«("+^^^ J2 f?^""^"''^ {imn - 2h^) A^{i^ + m7r/26). (6.4) 

meZ 

In particular, if q2{x) is a constant on the interval [0, 6], say q2{x) — 7, from (6.2) we 

get 

Mk) = -^2{^-^-'''')^ Bsik) = ^2 - ^-"""') ■ (6-5) 
Conversely, if A2{k) and B^{k) are as in (6.5), then with the help of 

"OO Aka 



— — dk = i7rsgn(a), 

-00 



where sgn(a) denotes the signature function, we recover q2{x) — 7 via (6.3). Alternatively, 
q2{x) can be recovered by using (6.4). 

In another particular case, in which q2{x) = e"^ on the interval [0,6], where a is a 
constant, from (6.2) we get 

g(a— 4ife)6 2 g(a— 2ifc)6 ^ 

'^'^^^= 2ik{a-Aiky '''^^^ = - 2ik{a-2iky ^'"'^ 
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Conversely, if A^{k) and B^{k) are as in (6.6), then we recover q2{x) = e°-^ via (6.3) with 
the help of a contour integration. Alternatively, we can recover q2{x) by using (6.4). 
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